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Abstract 

The conflicting arguments given in the discussion forum of Metamaterials 20 11 on the possible forms 
of macroscopic Maxwell equations are lead to a convergence by noting the relationship among the 
employed material variables for each scheme. The three schemes by Chipouline et al. using (A) stan- 
dard P and M (Casimir form), (B) generalized electric polarization Pll (Landau-Lifshitz form), 
(C) generalized magnetic polarization Ma (Anapole form) are compared with (D) the present au- 
thor's scheme using standard current density J. From the reversible relations among the transverse 
components of these vectors, one can easily rewrite one scheme into another The scheme (D), the 
only one among the four providing the first-principles expressions of susceptibility and also leading 
to a non-phenomenological Casimir form in terms of the four generalized susceptibilities between 
{P, M} and {E, B}, is concluded to be a more natural form than (B) and (C) as a single suscepti- 
bility theory. 

1. Introduction 

In the conventional macroscopic Maxwell equations (M-eqs), the variables of matter are usually rep- 
resented by the electric and magnetic polarizations P and M, in contrast to the microscopic M-eqs 
where we need only current (and charge) density J (p). In view of the general relationship J = 
dP/dt + V X M, the description in terms of P and M uses redundant variables. The constitutive 
equation to be required should be a single equation relating a single vector quantity of matter with that 
of EM field (E or B, or else). This aspect has been known for a long time, but not satisfactorily worked 
out. Landau and Lifshitz (LL) HI proposed to use the M-eqs V x E = —dB/dt, V x B = dD/dt, 
(i.e., H = B). This means that one uses a new variable Pll, defined by J = dPLi/dt, containing 
the both characters of P and M [2, 3|. The constitutive equation in this case relates D and E through a 
single susceptibility. LL discuss its symmetry properties, but not its quantum mechanical expression. 

Recently, Chipouline et al. (CST) f4l considered the possible forms of macroscopic M-eqs obtained by 
macroscopic averaging of microscopic M-eqs. Noting the non-uniqueness in the equations V • P = —p 
and V • J = —dp/dt, i.e., the fact that P and J may contain Vx of an arbitrary vector function, they 
consider three choices of matter variables, (A) usual P and M (Casimir form), (B) Pll (LL form), 
(C) Ma defined by J = V x Ma (Anapole form). They discuss the relationship among the three 
cases, including the possibility of their mutual transformation, but no quantum mechanical consideration 
is given about the form of susceptibilities. 



Another single susceptibility theory for macroscopic description was developed by the present author in 
terms of current density J tSJ. This theory leads to the constitutive equation 

J{k,uj) = Xcm{k,uj) ■ [A{k,uj) + {c/iuj)EextL{k,uj)] (1) 

where A is the transverse (T) vector potential in Coulomb gauge, and .EcxtL the longitudinal (L) electric 
field due to external charge density. (The L field due to internal charge density is taken into account as 
the Coulomb potential in matter Hamiltonian.) The macroscopic susceptibility Xcm is derived via long 
wavelength approximation of the microscopic (nonlocal) one as 

Xem{k,u}) = V'^[g,y{uj)iou{k)iuo{-k) + K{u)ii,o{k)io,y{-k)] , (2) 

V 

where v (and below) is the quantum number of matter eigenstates, V the quantization volume of k, 

Quiui) = 1/ {Ej^o — huj — iO+ ) — 1/ Ei^o ,hu{uj) = 1/ {E^q + fiuj + iO+ ) — 1/ EyQ , and E^q the excitation 



energy from the ground state. The matrix element of current density can be written as 

V(fc) = (exp(-ifc • r)/V)[J^, - ik ■ qJ=,') +ikx M^, + 0{k^)] (3) 

J,, = y dr {f,\JoW), M,, = Mjr^ + M^;i^^ , (4) 

k ■ Q^lly^ = XI / {< A'K'^'? - • Pi> ^('^^ -r) + 5{ri - r) {ri - f)k ■ pg\u >} , (5) 



where Joir) = J2e{^e/'^''^e){Pi ^(.fe — ''')+ 5(r^ — r)pi}, and L^^r) = (r^ — r) x p^, is the angular 
momentum of the ^-th particle with respect to the center coordinate f of the (/i, v) transition to make 
Taylor expansion of i^y{k). The zero-th and first order moments J ^j^y, Q^TJ^ , are nonzero for elec- 
tric dipole, electric quadrupole and magnetic dipole transitions, respectively. This result covers all the 
cases of linear response, including chiral susceptibility (Hill. This scheme should be added to the list of 
CST as the fourth item (D) specified by the use of matter variable J. It may be called "natural form", 
since it does not use unfamiliar variables as in (B) and (C). 

A discussion forum was held in the Metamaterials 2011 (Barcelona) about this problem. The large 
number of participants shows the general interests in this very fundamental problem. The discussions 
were rather conflicting with premature arguments, since the participants had not been well informed 
beforehand about the contents of other parties. Later the present author made a visit to have more 
detailed discussions with the CST group, which has resulted in this article unifying the schemes (A, B, 
C) with (D). 

2. Unification of the four forms 

The essential point for the unification of different schemes is to note that the non-uniqueness introduced 
by V • P = —p and V • J = —dp/dt is only for the T components of P and M, which is because 
the two equations give constraint only to the L components. This requires a refinement in the defining 
equations of the CST's classification, i.e, instead of J = dPiL/dt = V x Ma we should use 

j(^) = dP^^l/dt = V X M^'p . (7) 

This means that the L component of J is common to all the schemes (A, B, C, D). Namely, the choice 
of the variables should be 



(A) Casimirform: pC^), M^^) and j(^) , 

(B) LL form : and J^^^ , 

(C) Anapole forni : M^p and J^^^ , 

(D) Natural form : jC^) and j(^) . 

In (A) and (B), J^^-* may be replaced by p(^) and P^^\ respectively, which are equivalent to (i/tj) J^^\ 

Equation (|7]) can be solved as 

P^j^j^{k,uj) = {i/uj) J^'^\k,uj), MP{k,uj) = {i/k'^)k X J^^\k,uj). (8) 

This means that the constitutive equation for J, already known in the scheme (D), can be transformed 
into those for P^^ and Ma as 

pS2 = {i/u;)xcm-[A-ii/uj)E,^tL], (9) 
= {i/k^)k X [Xcm • {A - {i/uj)E,,tL}F^ . (10) 

The transformation from (D) to (A) is discussed in ||5] (Chap. 3) and ||6l by using the explicit expression 
of Xem> which defines the four susceptibilities, {(xeEi Xob)) (XmEi XmB)}> i-e., the electric and magnetic 
susceptibilities induced by E and B. This rewriting is reversible if one uses the microscopic expression 
of 

Xcm- Namely, the four susceptibilities can be put together to form the single susceptibility Xcm- In this 
sense, "the Casimir form derived from (D)", to be called scheme (A)*, is a single susceptibility theory. 
But the Casimir form with phenomenologically determined susceptibility, scheme (A), has no guarantee 
to be a single susceptibility theory. 

The argument given above shows that the four schemes can be transformed to one another From an 
arbitrary constituive equation one can derive all the other ones. This means that the dispersion equation 
should be same for all the schemes, i.e., the one already known in (D) 

det\[e-{-f]l-f,oX^^{k,oj)\=0 (11) 
c 

can be used also for (A)*, (B), and (C). 

To sum up, the scheme (D) is conceptually the simplest and practically the most informative scheme at 
present among the possible forms of macroscopic M-eqs. 

The author acknowledges the discussions with Dr A. Chipouline and S. Tretyakov for improving his 
understandings about the schemes (A), (B), and (C). 

References 

[1] L. D. Landau and E. M. Lifshitz, Electrodynamics of Continuous Media, 1st ed.(Pergamon Press, Oxford, 

1960), Sec.83, (Sec. 104 of 2nd ed., 1984) 
[2] Yu. A. H'inskii and L. V. Keldysh, Electromagnetic Response of Material Media (Plenum Press, New York, 

1994), Sec. 1.3 

[3] V. M. Agranovich and V. L. Ginzburg, Crystal optics with Spatial Dispersion, and Excitons, (Springer Verlag, 
Berlin, 1984), Sec.6, 

[4] A. Chipouline, K. Simovski, and S. Tretyakov, Metamaterials 6 (2012) 77-120 

[5] K. Cho, Reconstruction of Macroscopic Maxwell Equations, Springer, Heidelberg, 2010 

[6] K. Cho Proc. Metamaterials 201 1, Barcelona, p.672-674 



